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Abstract 

Aim of this work was to test the idea of J.-Y.Ollitrault about the 
parallel squeeze-out type of transversal flow in mid-rapidity region Q. 
For this purpose we have performed a computer simulation of the ex- 
panding pion gas created in non-central Pb-Pb 158GeV/n collisions. 
A squeeze-out type of asymmetry parallel to impact parameter in az- 
| imuthal distribution of pions is found and studied. The asymmetry is 

explained as a consequence of geometry of non-central collisions and the 
^ rescattering process. 

o : 

\0 ■ PACS numbers: 25.75 

. 

■ , 1. Introduction 

Ph. 

CD . 

^ ■ Azimuthal asymmetries in transverse momentum distributions of particles 

measured in relativistic heavy ion collisions (HIC) were observed for the first 
time by plastic-ball detector in Berkeley ||. Since that pioneering experiment 
the asymmetries in relativistic HIC were measured also at higher energies [^|, ffl . 
Typical types of asymmetries (bounce-off and squeeze-out) are explained as a 
consequence of collective behaviour of nuclear matter. This understanding of 
the origin of azimuthal asymmetries in non-central HIC succesfully explains 
most of the experimental data. 

However at AGS experiment with Au+Au 11.4 GeV/n collisions E877 col- 
laboration reported about the unclear origin of a new type squeeze-out ef- 
fect parallel to impact parameter. It seems that E877 collaboration has found 
the first experimental indications of the effect predicted by J.-Y.Ollitrault [ffl. 
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This type of squeeze-out effect results from the interaction of produced parti- 
cles among themselves. 

In this work we investigate asymmetry in transverse momentum distribu- 
tion of pions created in the process of HIC. No collective behaviour of nuclear 
matter is supposed in the calculation. Asymmetry results from the geometry 
of non-central HIC and the rescattering process. 

Paper is organized as follows: In Section 2 we give an intuitive explanation 
of the origin of the asymmetry. Description of the computer simulation is given 
in Section 3. In Section 4 we present results of our simulation - dependence the 
asymmetry on impact parameter and also other features of the studied effect. 
We finish this work with short summary and conclusions. 

2. Asymmetry as a consequence of rescattering 

Usual explanation of the azimuthal asymmetries in momentum distribu- 
tions of particles in HIC is based on the collective behaviour of nuclear matter. 
Hydrodynamical models 0, || describe the collision process of HIC using the 
equation of state for nuclear matter. Substantial assumption of these models is 
a thermalisation process. As we shall see the existence of thermal equilibrium 
is not required for the existence of our type of asymmetry. 

Number of secondary particles created in HIC at present SPS energies 
(158GeV/n) is substantially larger than the number of nucleons contained in 
the colliding nuclei. Therefore the situation is different than that at lower 
energies and the interaction of the secondaries among themselves becomes 
significant. As a signature of the rescattering process among the produced 
pions we regard the pt dependence of transverse size parameters Rt [0] extracted 
by HBT technique ||. The decrease of R t with increasing p t was explained 
as a consequence of the rescattering of pions in the simulation [[| where the 
central S-Pb 200 GeV/n collisions were studied. 

We think that the rescattering of secondaries can demonstrate itself in the 
experimental data also as a transverse flow of pions in non-central collisions. 
For the existence of this effect no collective behaviour of nuclear matter is 
necessary. Therefore the phenomenon can be observed also in the future HIC 
experiments (RHIC, LHC) where nuclear transparency region is expected to 
be reached. 

In this section we explain the origin of the predicted transverse flow of 
pions as a consequence of the geometry of non-central HIC and the rescatter- 
ing process. On Fig. la we show the geometry of non-central A+A collision in 
transverse plane. Secondary particles i.e. also the pions are created mainly in 
the overlapping region where nucleon-nucleon collisions happen. Shape of the 
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overlapping region depends on impact parameter b and is azimuthally asym- 
metric in transverse plane. Since we do not consider any collective behaviour 
of nuclear matter it is natural to assume that the initial distribution of pions 
in transverse momenta ^(pt) is azimuthally symmetrical: ^ s (p t ) = ^(\pt\). 
Because of the simplicity of explanation we do not write longitudinal compo- 
nents of p and x in this section however the simulation described in the next 
section is performed in 3 dimensions. 

Denoting the space distribution of the points of creation of pions by &(x t ) 
we have the following initial condition for the pion gas created in HIC: 



X{x t ,p t ) = $ A {x t )-* s (p t ) 



(1) 



Because of the rescattering process the original x — p non-correlated distri- 
bution (1) becomes x — p correlated [HJ and (as we shall see) the asymmetry 
in $ (x t ) leaks into asymmetry in the resulting transverse momentum distri- 
bution. 

^ s {p t )-<$> A {x t )->X A {p u x) (2) 

where A denotes asymmetry in distribution X(p t ,x) in momentum. 

This effect can be understood in the following way: Let us have two groups 
of pions: A group of parallel pions with the momentum parallel to b and the 
group of orthogonal pions with the momentum orthogonal to b. Because of the 
asymmetry in Q A (xt) the probability of collision of the orthogonal pion is (in 
average) higher than the probability of collision for paralell pion (see Fig. lb). 
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Fig.l Geometry of non-central collision in transversal plane 

The reason for this is very simple: Number of pions located in the direction 
orthogonal to b is higher because of the bigger size of initial "cloud" of pions 
in this direction. Thus during the rescattering process orthogonal pions collide 
more likely than the parallel pions and consequently collisions of type Fig.lc 
are more frequent than the collisions of type Fig. Id. Such conditions lead to 
the excess in the number of paralell pions. 
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This is undoubtedly a non-equilibrium process. Resulting asymmetry in 
the momentum distribution of pions freezes in the considered pion gas because 
of the expansion of the system. 

Azimuthal asymmetry in transverse momentum distribution of pions after 
the rescattering process can be measured by detectors as the excess of pions 
in the direction paralell to the impact parameter. As we shall see in the 
next sections predicted asymmetry should demonstrate itself as a second order 
asymmetry in the fourier analysis of transverse flow 
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On Fig. 2 we show the result of a toy simulation of the effect described. 
Initial distribution of pions in momenta was taken from the central S-Pb 200 
GeV/n HIC simulation [10, [| and therefore it was azimuthally symmetric. 



Initial distribution of pions in x t space is artificially asymmetric as it is shown 
on Fig. 2b. After 20fm/c of the time evolution the rescattering process leads 
to the asymmetry in transverse momentum distribution as it is clearly seen 
from Fig. 2c. 
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Fig. 2 Toy simulation of the effect, a) is pt distribution before the rescat- 
tering. c) is pt distribution after 20fm/c of rescattering process. Cut 
in transverse momentum pt > 300MeV is applied to enhance the visual 
effect of the asymmetry, b) is xt distribution before the rescattering 
process and d) after the 20 fm/c of expansion. 

Simulation of the expanding pion gas for "real" Pb-Pb 158 GeV/n non- 
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central collisions is described in the next section. 



3. The Simulation 

For the simulation of the expanding pion gas created by "real" Pb-Pb 158 
GeV/n non-central collisions two independent programs were used: 

1) Cascading generator [O] which generates the initial momenta and posi- 
tions of pions as a result of independent nucleon-nucleon collisions. 

2) Rescattering program JTO] which simulates time evolution of the inter- 
acting of pions created in HIC. 

First we shall describe main interface structure of the simulation shown 
on Fig.3. Pb-Pb 158 GeV/n collisions were simulated by the cascading gen- 
erator (CG) for random orientations and selected values of impact parameter 
b. Information about initial momenta, time and place of the creation of pions 



tt(x,p, t) produced by CG was used as input for the rescattering program JlO 



Final momenta of the interacting pions were selected from the output of the 
rescattering program and analyzed for transversal asymmetries. Orientation 
of the impact parameter was not supplied by CG and its determination is 
described in subsection 3.3. 
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Fig.3 Main interface structure of the simulation 



3.2 Rescattering process 

Description of internal structure of the program used for the simulation of 



rescattering process can be found in [10]. This program was build according 



to the description of the simulation of central S-Pb 200 GeV/n collisions 



Results obtained by program |10[ for the central S-Pb collisions are in agree- 
ment with the results obtained in |§. We shall sketch here just main features 
of our program. 

Pions are treated as point-like objects in the simulation, position and mo- 
mentum of each pion is known during the simulation. Pions move in small time 
steps At = 0.1 fm as a free particles. Collisions happen if two pions appear 
to be at a distance smaller than the critical distance which is determined 
from the isospin averaged total elastic cross-section 
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D, 



(3) 



Momenta of pions after the collision are determined in CMS of the pair 
according to the isospin averaged differential cross section 



da/dQ = a(s) + b(s) cos 2 (6>) 



(4) 
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(Numerical values of the functions a(s), b(s) are calculated from data 

Then new momenta are transformed back to the global frame of the sim- 
ulation. Test for the relative distance is performed for every pair of pions in 
each time step. Not all of pions evolve in time at the beginning of simulation. 
Cascading generator produces also the information about the time of creation 
for each pion therefore pions are tested for the existence in the global frame of 
simulation. 




Fig. 4 Dynamics of tttt collisions. Collisions can happen only if both 
pions are allowed to interact (solid lines). 

Moreover existing pions are restricted from the interaction for time Tj - 
formation time [15] after their creation and also for time after each collision 
(see description in []T0|] ) . These two phenomeno logical parameters allow to 
influence total number of collisions in the simulation and also some features of 
the expanding of pion gas (see Tab.l). 



3.3 Data analysis 

Events of the exact values of impact parameter b = 3,5, 7, 9, 11/m were 
used for the analysis. Approximate orientation of impact parameter was de- 



termined from output of the cascading generator. Cascading generator ||T2 
simulates A-A collision with the target nucleus at the center of coordinate 
system. Therefore the initial positions of pions which lie mainly in the over- 
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lapping region of nuclei are shifted from the center of coordinate system in the 
direction of impact parameter (see Fig. 5). 

A procedure close to Danielewicz - Odyniec method |16|] was used for deter- 
mination of the approximate orientation of impact parameter. For each event 
vector 



<3 = E 



(5) 



{x\ are initial positions of created pions in transversal plane) was constructed 
and its orientation was used as the orientation of impact parameter b. Final 
momenta of pions after the rescattering process were rotated to have the same 
orientation of the approximate impact paramenter. 



10 

5 

-5 
-10 



: ; oo*oo- 



-10 -5 5 

b = 3 fm ROT 



10 



10 
5 

-5 
-10 



-10 



-5 5 

b = 7 fm ROT 



10 



10 
5 

-5 
-JO 



-10 -5 5 10 

b = 7 fm NoR 



Fig. 5 Initial positions of pions in transversal plane obtained from the 
cascading generator. 10 events for b=3fm and b=7fm are rotated to 
have the same orientation of b. On figure b=7fm NoR the rotation is 
not performed. 

Rotated momenta of pions were used for fourier type of analysis of transver- 
sal flow 
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For each pion azimuthal angle of momentum in respect to the 
impact parameter was determined and added to histogram of azimuthal dis- 
tribution R((f>) (see Fig. 6). Because of the known character of asymmetry and 
the rotation of events into the direction with the same orientation of impact 
parameter we have fitted the normalized histograms R ((f)) to the function 



R n ((j)) = 1 + S 2 cos(20) 
Results of the analysis are presented and discussed in the next sections. 

5. Results 



(6) 



First the presence of asymmetry was tested on the set of 24 artificial S- 
Pb 200 GeV/n asymmetrical events. A strong final asymmetry in azimuthal 
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distribution of pions in transversal momentum is visible directly from Fig. 2. 
We have tested also dependence of the strength of the effect on the dynamics of 
rescattering process. Normalized histograms of R(4>) distribution were fitted to 
function (6) for different values of parameters Tf and Tj. The results obtained 
are summarized in Table 1. 



T f /Ti[fm] 


0.2/0.2 


0.2/0.5 


0.2/0.8 


0.5/0.2 


0.5/0.5 


0.8/0.0 


1.0/0.0 


S 2 ± 5% 


0.43 


0.36 


0.42 


0.42 


0.41 


0.36 


0.29 


N C oii/7t 


4.4 


2.6 


2.1 


1.9 


1.3 


0.9 


0.5 



Tab.l Results of fit for 24 artificial S-Pb events. ^coii/' K is average 
number of collisions per tt during the rescattering process. Total number 
of pions in one event was 700. 

Then "real" non-central Pb-Pb 158 GeV/n collisions were studied using the 



output of cascading generator ||l2"f . We have run the program for five sets of 10 
events with impact parameter b = 3,5, 7, 9, 11 fm in order to study dependence 
of the asymmetry on the impact parameter and additional 20 events for b = 
7fm were run in order to verify whether the statistics of 10 events per set 
was sufficient^. Pions in rapidity range (—1, 1) were selected for the analysis. 
Events were rotated to have the same orientation of impact parameter using 
the procedure described in section 3.3. On Fig.6 we show the histograms of 
distribution R((p) for the ROT - rotated (to the same orientation of b) and 
non-rotated (NoR) events for b=3fm and b=9fm. Normalized distributions 
R N (4>) were used for the fit. 
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2 For b — llfm 20 events were run because of low multiplicity of pions in these events. 
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Fig. 6 Histograms of the azimuthal distributions of pions for b=3fm and 
b=9fm. Numbers of particles in the bins of histogram for b=3fm are 
higher than in the case b=9fm because of the higher total multiplicity of 
pions in the collisions b=3fm. Normalized histogram R N (4>) for b=3fm 
is artificially shifted up. 

Results of the fit are summarized in Tab. 2. Simulation of 20 additional 
events for b=7fm showed that the statistics of 10 events is sufficient for the 
qualitative analysis we have performed. Main result of our simulation is visible 
directly from Fig. 7: Asymmetry increases with impact parameter in range 
3 — 9fm and it falls down for b > 9fm. Turnover point of S2 lies between 7 
and 11 fm. We think that the position of the turnover point depends on Tf, 
parameters used in the simulation. 

Rapidity dependence of the asymmetry was analyzed for events with impact 
parameter h—lfm. Our small statistics allowed us to perform only a very rough 
analysis of rapidity dependence of S2 coefficient. Results in the following table 
are averaged for forward and backward rapidities arround central rapidity what 
is in our simulation Y c = 0.0. 





Y 




(0,1) 


(1,2) 


(2,3) 


s 2 


0.14 


0.15 


0.9 



We have tried to find also first order asymmetry (bounce-off) in our data 
which would demonstrate itself as a non-zero value of Si coefficient in the 
fourier expansion: 

R{<j>) = A ■ (1 + S 1 ■ cos(0) + S 2 • cos(20)) (8) 
For events b = 7fm we have obtained following results: 



Y 


(-3,-1) 


(1,3) 


Si before 


+0.007 ±0.01 


-0.003 ±0.01 


Si after 


+0.016 ±0.01 


-0.02 ±0.01 


(p x ) in MeV 


2.22 


-1.96 


(p y ) in MeV 


0.06 


-0.06 



Si before is calculated from the output of cascading generator and Si after 
was computed for events after the rescattering process (see Fig. 3). Slight 
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signature of the first order asymmetry is present in the output of cascading 
generator what can be a consequence of the cascading of nucleons in spectator 
matter included in the generator (T^|. However this effect is not above the 
statistical errors in our set of evets h—lfm. 

Calculation of average (p x ) and (p y ) of pions in the resulting events b=7/m 
confirms that the first order asymmetry is very weak in our data. 

We have computed R p parameter used in the study of squeeze-out effect at 



BEVALAC/SIS energies |TJ 



P ~ (Pi) - (Py) 2 { ) 

where x direction is parallel to impact parameter. For b=7/m our R p of rotated 
(ROT) and non-rotated (NoR) events is: 

R* OT = 1.39 R» oR = 0.97 (10) 

An interesting behaviour of asymmetry coefficient S2 was found in our data 
analysis. As it is shown in Tab. 2 coefficient S2 is significantly higher for high 
p t pions (p t > 300MeV). 



b 


3fm 


5fm 


7 fm /30 


9fm 


11 fm/20 


(P, > 0) S 2 


0.06 ±0.01 


0.09 + 0.01 


0.135 ±0.010 


0.20 ± 0.02 


0.13 ±0.02 


(P t > 300) S 2 


0.13 


0.15 


0.21 


0.30 


0.21 


N co ii/7r 


4.5 


4.2 


3.4 


2.6 


1.1 



Table 2. 



Tab. 2 Asymmetry coefficients for Pb-Pb 158 GeV/n non-central events. 
7 fm/30 means that the fit was performed for 30 events. 

We have tested whether this behaviour is present also in E t sensitive anal- 
ysis. For this purpose we have filled histograms of R(4>) distribution by the 
value E t for each pion. P t dependence of S2 coefficient was still present in the 
results. 



5. Summary and Conclusions 



We have studied asymmetry in the azimuthal transverse momentum distri- 
bution of pions for non-central Pb-Pb 158 GeV/n collisions. For the presence 
of this type of asymmetry no collective behaviour of nuclear matter is nec- 
essary. It is a consequence of the geometry of non-central collisions and the 
rescattering among the pions. 
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The simulation showed that the asymmetry increases with the impact pa- 
rameter in the range up to 9fm. Turnover point is located between 7 — 11 fm 
in our data. Decrease of asymmetry coefficient S2 for b = 11/m is most likely 
a consequence of low number of pions participating in the rescattering process 
combined with small total size of the overlapping region in comparison to for- 
mation path Lf ~ c-Tf of pions. Non-interacting pions in the "formation time 
stage" do not feel the asymmetric shape of the overlapping region of colliding 
ions and therefore the asymmetry decreases. If this is the main reason of the 
decrease of asymmetry at peripheral collisions then the position of the turnover 
point can be sensitive to the value of formation time parameter Tf used in the 
simulation. 

We think the effect studied in this work was already confirmed experimen- 



tally by NA49 collaboration [18| using Ring Calorimeter setup [[19 . 

Unexpected p t dependence of the S2 asymmetry coefficient was found in 
the analysis of our data. This seems to be an interesting prediction however 
the origin of this pt dependence is not clear at present. It can be a consequence 
of our scenario of rescattering process and therefore it does not need to appear 
in experimental data. Results obtained from NA49 TPC could answer this 
question. We hope some data will be presented during QM'96 conference. 
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Fig. 7 Dependence of S2 coefficient on impact parameter. For high-pt 
pions the 52 coefficient is systematically higher than in the case of full 
Pt range pions. 

At the future HIC experiments the multiplicities of secondaries will be much 
higher. Because of the nuclear transparency phenomenon the mechanism re- 
sponsible for the transversal flow of nucleons at lower energies (up to lOGeV/n) 
can play a little role. In this case the studied effect can be substantial for the 
transversal flow phenomenon at HIC on RHIC and LHC. 

Acknowledgements 

Author is indebted to prof. J.Pisiit for supervising this work. Special 
thanks are directed to P.Zavada - the author of the Cascading Generator for 
the possibility to use output of his program. We express our gratitude to 
member of NA49 collaboration T.Wienold for the important comments and 
suggestions. 

This work was originally motivated by the discussion with J.-Y.Ollitrault. 
We thank to J.-Y.Ollitrault also for the subsequent fruitful discussions. 



12 



References 

J.-Y.Ollitrault, Phys.Rev. D 48 1132 (1993) 

H.A.Gustafsson et AL., Phys.Lett. B 142 141 (1984) 

J.Barrette et al., Phys.Rev. Lett. 73 2532 (1994) 

E877 Collaboration in Proceedings of Quark Matter '95 

J.D.Bjorken, Phys.Rev. D 27 140 (1983) 

J.-Y.Ollitrault, Phys. Rev. D 46 229 (1992) 

NA 44 Collaboration Phys. Rev. Lett. 74 17; 3340 (1995) 

D.H.Boal et AL., Rev. of Mod. Phys. 62 553 (1990) 

T.J. Humanic, Phys.Rev. C 50 2525 (1994) 

P.Filip, Acta Physics Slovaca Vol 50 9 (1996) 

S.Voloshin, Y.Zhang in hep/ph/9407282 

P.Zavada, Phys.Rev. C 42 1104 (1990) 

M.Prakash et al., Phys. Rep. 227 321 (1993) 

P.Estabrooks, Nucl.Phys. B 95 322 (1995) 

J.PlSUT et. AL. Zeit.f.Phys. C 67, 467 (1995) 

P.Danielewicz, G.Odyniec Phys.Lett. 157B 146 (1985) 

H.H.Gutbrod et AL., Phys.Rev. C 42 640 (1990) 

T.Wienold et AL., LBL 38506, to appear in Proceedings of 12th Winter 
Workshop on Nuclear Dynamics, Snowbird 1996, Utah, USA 

[19] S.Margetis et al., Nucl.Phys. A 590 355 (1995) 



13 



